absTol = eps()^(2/3);

x = tseries(qq(1,1):qq(3,2),sin(1:10)) ;

observ = double(hpf(x)) ; 
expect = [ ...
    0.332158721583314
    0.287235347551190
    0.242630293683581
    0.199050668944550
    0.157140138369643
    0.116944958766749
    0.077813846685613
    0.038847793390374
   -0.000490226884807
   -0.040143170872499 ...
   ] ;
myassert(observ, expect, absTol) ;

observ = double(hpf2(x)) ;
expect = [ ... 
    0.509312263224583
    0.622062079274492
   -0.101510285623714
   -0.955853164252478
   -1.116064413032781
   -0.396360456965675
    0.579172752033176
    0.950510453233008
    0.412608712126564
   -0.503877940016871] ;
myassert(observ, expect, absTol) ;

